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Abstract 

We present results of our studies of various scattering properties of topological 
solitons on obstructions in the form of holes and barriers in 1+1 dimensions. Our 
results are based on two models involving a ip 4 potential. The obstructions are 
characterised by a potential parameter, A which has a non-zero value in a certain 
region of space and zero elsewhere. In the first model the potential parameter is 
included in the potential and in the second model the potential parameter is in- 
cluded in the metric. Our results are based on numerical simulations and analytical 
considerations. 



1 Introduction 

Some surprising aspects of the scattering of solitons have been attracting the attention 
of researchers and have been studied in many papers [1-3]. These studies have shown 
that a topological soliton behaves like a classical point particle when it is scattered on 
a potential barrier. In such a case the soliton follows a well denned trajectory and as it 
meets the barrier it slows down and, if it has enough energy to 'climb' over the barrier, it 
gets transmitted otherwise it gets reflected with almost no loss of energy. Hence, soliton's 
scattering on a potential barrier is very elastic and thus the soliton behaves like a classical 
point particle. However, a soliton displays more interesting behaviour when it encounters 
a potential hole. Consider a classical point particle encountering a potential hole. Then, 
when the particle enters the hole, the hole speeds it up, no energy is lost, and the point 
particle always gets transmitted. However, at low velocities, below a critical value, the 
scattering of solitons on a potential hole exhibits some aspects of a quantum like behaviour 
in the sense that, sometimes, the solitons get reflected by the hole although most of the 
time they get trapped in the hole. Of course such reflections are not possible for a classical 
point particle. A soliton with enough energy, ie with velocity above a critical value can 
get out of the hole but it comes out with much reduced velocity as during the whole 
process it radiates some of its energy [1-2]. 
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on an obstruction in two Lagrangian models. The Lagrangians are so constructed that, 
far away from the obstruction, the soliton in the two models looks exactly the same. 

In the numerical part of this work we used square well potential barriers and holes of 
width 10. The simulations were performed using the 4th order Runge Kutta method of 
simulating the time evolution. We used 1201 points with the lattice spacing of dx = 0.01. 
Hence, the lattice extended from -60 to 60 in the ^-direction. The time step was chosen 
to be dt = 0.0025. We used the absorbing boundary conditions. 

1.1 Model 1: X(x) ip 4 potential. 

First we consider the model, in (1+1) dimensions fx—0,1, defined by 

t=±d lt <p& l <p-U(<p) (1) 

with the potential 

U(<p) = \(x) (y? 2 -l) 2 , 

where A = Ao + A(x), and X(x) is a potential parameter which has been inserted into the 
Lagrangian to take into account the effects of obstructions, holes and barriers, and so is 
nonzero only in a certain region of space. 

After varying the action, the equation of motion is 

dpdPip + 4A (x) ip (if 2 - l) = 0. (2) 

This equation cannot be solved analytically because the potential has a spatial de- 
pendence. However, this equation can be solved numerically by placing the soliton in the 
region where A = by making and then evolving it from there. 

Thus the potential parameter \(x) has a constant value, Ao, in the region where the 
soliton is located and far away from the position of the soliton \(x)ip (y? 2 — 1) is zero as 
if ~ 0. Moreover, as the system is fully relativistic, we can get a time-dependent solution 
by simply boosting the static one. This gives us a soliton moving with velocity u. Thus 
(assuming that x is nonzero and t is small we can put 



(p (x, t) — ± tanh yyy 2X (x — x — ut) j , (3) 

where 



7 = -T^r (4) 



u 



Here the plus and minus signs represent the kink and anti-kink solutions respectively. 
At t — 0.0 the soliton is far away from the barrier/hole and this description described the 
soliton very well. In what follows we put Ao = 1.0 and so our initial conditions for the 
field and its time derivative are given by 

(p(x,0) = ±tanh(7v / 2(x-Xo)) (5) 
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o ip (x, 0) = =F wyv 2 seen (jV2{x — xo)J (6) 
The energy density of the field is given by 

t=\{d^f + \{d lV f + ( V 2 -if . (7) 

1.2 Model 2: A \(p A model with a position dependent metric. 

In this paper, we consider also the behaviour of a soliton in (1+1) dimensions when the 
potential parameter (A(x)) arises in the space-time metric, ie the metric is given by 



1 + A (x) 
-1 



Thus we are putting the obstruction into the model in the same way as was done in 
[3] in the sine-Gordon case. 
The action is now given by 



£ (ip, dpip) dx dt, (8 
where £ is the Lagrangian density and g is the determinant of the metric, ie g (x) MI/ 



(9) 



with the potential still being given by 

U{ V ) = (^-lf. 

Note that the coefficient of the potential in U(<p) is chosen to be 1, ie it is the same 
as in model 1 away from the obstruction. The equation of motion is now: 



[l + \{x))dl V -dlip 



-<9iA(x)<9i<^ + 4^( 



0. 



(10) 



2 | 1 + A (x) , 

The topological charge is independent of the metric and does not depend on the 
potential. The equation of motion (1101) cannot be solved analytically because of the 
spatial dependence of the potential parameter. But, like before, a solitonic solution can 
be found by making the potential parameter constant, A (x) = Ao, and then use that 
solution as an initial condition to solve the equation (TTUT) numerically. Note that the 
equation (flOil can be rewritten as 



12) 



(1 + A ) d V - %tp + 4 (1 + A ) if {if 2 - l) = 0. 
And the approximately time-dependent solution for the field is given by 



ip (x, t) — ± tanh 



vv 1- u 2 (l + X 0) 
At t = 0.0, Ao = 0.0 and so the initial conditions are 



X — Xq — Ut) 
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Figure 1: The potential barrier as seen by the soliton in the two models 



ip (x, 0) = ± tanh 



1 _ U 2 (1 + Aq ) 



(x -x )\, 



d O (p(x,0) = TU\ 



sech 



, 1- u 2 (l + X ) V V 1 - ^ 2 (l + Ao) 
which are, in fact, the same as the initial conditions of model 1. 



(x - x ) , 



(13) 
(14) 



The energy density is then given by 

7 00 
2 



a 00 1 / x ' 

^(^) 2 + -(^) 2 +(^ 2 -l)' 



(15) 



2 Barrier Scatterings in both models. 



The scattering of a topological soliton on a potential barrier, in model 1, was found to 
be nearly elastic. However, in model 2, the scattering was found be less elastic as, in 
this case, the soliton radiated away a small amount of its energy, although this radiated 
energy was very small (almost negligible) when the velocity of the soliton was far below 
its critical value. However, it increased as we sent solitons with higher velocities. 

The shape of the potential barrier, as seen by the soliton, can be found by placing the 
soliton with a zero velocity at different positions and calculating its total energy. This 
has given us figure 1 which shows how a soliton would see a potential barrier of height 
0.5 in both models 

However, performing numerical simulations we have found a difference between the 
two models. In model 1 the behaviour of a topological soliton resembles a classical point 
particle as only a very small amount of energy is being radiated out during the scattering 
process. However, in model 2, this amount of radiation is much larger and it increases as 
we increase the height of the barrier. Figure 2a and 2b shows the time evolution of the 
position of the soliton (moving with velocity 0.8) for different barriers heights in models 1 
and model 2 respectively. We see very clearly that in model 2 the soliton loses more energy 
(through radiation) and consequently leave the scattrering region with lower velocity. Not 
surprisingly, this effect grows with the increase in barrier's height in both models but is 
much more pronounced in model 2. 

We have also found a small variation of outgoing velocity of the soliton as a function 
of its incoming velocity; again this effect is more pronounced in model 2. Figure 3 shows 
our results for a potential barrier of height = 0.5. Figure 4 indicates, how in model, this 
variation Oscillation) get more pronounced with the increase the height of the barrier. 
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Model 1 : position of soliton vs lime at different potential heights 



Model 2 : position of soliton vs time at different potential heights 
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gure 2: The position of soliton vs time for various barrier's heights 



the outgoing velocity vs the incoming velocity for 
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gure 3: The v out vs Vi n for a barrier height of 0.5 in the two models 
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Figure 4: Model 2: v out vs v in for various barrier's heights 
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barrier's height of 0.5 



Figure 5: Ratio of the slopes of soliton's solutions vs incoming velocity in the two models 



2.1 Some Analysis of Our Results. 

We have found that the solitons behave in a similar way in the two models but that 
the differences from the behaviour of a point particle are more pronounced in model 2. 
This difference was found to be related to the oscillations of the kink corresponding to 
the change of its slope. These oscillations are generated because the slope of soliton 
changes as the soliton moves up a barrier. This encounter is somewhat sudden, the slope 
changes too much and so its starts oscillating. The soliton thus tries to readjust itself 
to the new value of the slope (corresponding to the new value of lambda) and during 
this readjustment some of its kinetic energy is converted into oscillating energy. To test 
such an interpretation we divided the slope of the soliton before the interaction with the 
barrier over the slope of the soliton when it is at the barrier and then have plotted this 
ratio as a function of the incoming velocity. Figure 5 shows the plot of the ratio for a 
barrier of height = 0.5 in the two models. From the figure on can see that below the 
critical value, it when the soliton gets reflected, the ratio of the two slopes is small and 
increases steadily as the soliton velocity increases. However, above the critical value, ie 
the soliton gets transmitted, the ratio in model 2 increases steadily but only slightly this 
time while in the model 1 it is almost constant. 

Furthermore, we have also studied the range of the oscillations and the average square 
of their frequency. We have found that above the critical value the average square of 
the frequency is directly proportional to the incoming velocity. This demonstrates that 
the energy radiated, which is proportional to the square of the frequency, increases as we 
increase the incoming velocity. This supports our argument that the soliton oscillates and 
that a fraction of its energy is converted into vibrating energy. 

The range of the oscillation is given by 

Sl ~ S2 (16) 



Siy S 2 — 1 

where Si is the slope of the soliton before and after the barrier. S 2 is the slope at the 
barrier. This result holds in the limit when the difference between the slopes is very small 
(so that we can make the approximation that sin(r5) ~ 5). The average of the frequency, 
u, is given by 

« = < 17 > 

where At is the time spent by the soliton at the barrier. Figure 6 illustrates the relation 
between the average square of the frequency and the incoming velocity, at heights 0.5 and 
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0.75. The figure shows that, as the incoming velocity increases, the oscillations increase 
and so does the radiated energy. 

We have found that the difference in oscillations and obviously in the behaviour of 
solitons in the two models is due to the different energy, ie mass, they have at the barrier. 
The more massive the soliton is the fewer oscillations it makes during the interaction. 
There is a greater sensitivity of the potential for less massive solitons than for the heavy 
ones. The mass of a static soliton, M rest) is given by the integral over the energy density: 



M, 



rest 



dx 



(18) 



3V% 

When the soliton is far from the obstruction A = 0.0 and so its energy ie mass 
M = 1.88557 in both models and this is exactly what was found numerically. For a 
moving soliton its is given by 



E 



J — < 



dx = 



M 



vT 



(19) 



where u is the speed of the soliton. 

Figure 7 shows a plot of the soliton rest mass as a function of the height of the barrier 
in the two models. 

We have also looked at the values of critical velocities in both models. In the table 1 
we give these values for various heights of the potential in both models. It is clear that 
the critical velocities in model 2 are greater than those of the model 1. That is very 
understandable since to reach the values of the rest masses in model 1, one needs a larger 
value of the gamma factor which thus leads to smaller critical velocities. 
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critical velocities 











AQ — V. 1 O 


Model 1 


~ 0.34 


~ 0.45 


~ 0.58 


~ 0.66 


Model 2 


~ 0.358 


~ 0.463 


~ 0.596 


~ 0.695 



Table 1: Soliton' s critical velocities for potentials of various heights 



:ls ci il vebcky vi Bsirior'i 'mic'i: ir two models 



Figure 8: Critical velocity as a function of barrier's height in the two models 



Figure 8 shows a plot of the critical velocity as a function of the height barriers in the 
two models. 

A soliton is an extended object. However, in a particle picture, the soliton rest mass 
energy is 

E= Mc 2 , with c= 1 (20) 

where M is the mass of the soliton. Thus, in the particle picture, if a soliton is sent 
towards a potential barrier with a critical velocity, u cr , then it would almost have enough 
energy to climb up the potential barrier. In this situation the soliton critical energy, E cr , 
has to be very close to the soliton rest mass. In a classical particle picture if the kinetic 
energy of the soliton is completely converted into the potential energy, the soliton will 
forever stay at the top of the barrier because the difference between the critical energy 
and the rest mass energy is zero. However, practically, the difference is very small and if 
we let the difference to be around ~ x where x is a very small number, then 



E r: 



M, 



rest 



X., 



where M rest is the rest mass of the soliton at the top of the barrier and we expect x to 
be very small. 

M , x 

E rr = , (21) 



uf. 



where M is the rest mass of soliton when it is far away from the potential which was 
found to be 1.88557. Thus, the critical velocity is given by 



u, 



M 



M res t + x 



(22) 



In model 2 the rest mass of the soliton at the barrier of height 0.5 was found to be 
2.3093, see fig 7. Putting x = 0.02 we estimate the critical velocity as 



1 ~ ( o.02+2 5 30 9 ) ~ - 592 ' 
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Figure 9: The potential hole as seen by the soliton 



which is in excellent agreement with what was found numerically (see fig. 8 where it is 
given by 0.596). 

3 Hole Scattering 

A classical particle is always transmitted through a potential hole. However, numerical 
simulations of the scattering of topological solitons [1-2] have shown that, when the po- 
tential hole is deep enough, there is a critical velocity below which the soliton gets trapped 
in the hole. Moreover, just below this critical velocity, the soliton can get reflected by the 
hole, thus showing behaviour which is more like of a quantum object. Of course, there 
is nothing 'quantum' about it - the soliton gets trapped in the hole and interacts with 
the radiation waves that it has sent out and ... occasionally, for very specific values of 
its velocity this interaction leads to the back escape of the soliton which lookes like its 
reflection by the hole. 

In the present work we have also performed many simulations of the scattering of 
solitons on the holes for the two models and have seen the same behaviour. Like in [1] 
we have found that the soliton get get trapped in the hole and that it can get reflected. 
In all cases it radiates some energy hence, if it comes out of the hole its velocity is lower 
than at its entry. 

In our work we placed a square hole of width 10 at the same position as the barrier in 
the previous simulations. Again, when the soliton is far away from the hole its rest mass 
M = 1.88557. When the soliton approaches the hole of -0.50 depth, it sees the hole as in 
figure 9 with no difference between the two models. In figure 10 we present a plot of the 
soliton rest mass as a function of the hole depth. 

As, when the soliton is in the hole, it has a lower rest mass in the model 2 than in 
the model 1, the soliton in model 2 has more 'spare' energy and so it it radiates it more. 
Figure 11 shows how the position of a soliton moving with a velocity of 0.65 evolves in 
time in the two models. From the figure we can see that the loss of energy in the model 
2 is higher than in the model 1. 

Incidentally the model 2 has one more interesting scattering property. We have found 
that for a hole depth just below -0.625, the solitons are always trapped. Hence there is 
no critical velocity. For example at a hole depth of -0.625 the critical velocity is ~ 0.97. 
Figure 12 shows a plot of the critical velocity as a function of depth of the hole for model 
2. 



9 



-0.6 -0.55 -0.5 -0.45 -OA -0.35 -0.3 -0.25 -0.2 
hole's depth 



Figure 10: Soliton rest mass vs hole's depth in the two models 




We have constructed two models involving two different ways of introducing a localised 
obstruction in the A<^ 4 model and on the scattering properties of solitons on obstructions 
in such models. In the first model the coefficient A was made to be a function of the space 
variable x, in the second one the obstruction was introduced through the modification of 
the metric. The scattering properties of solitons in both models were very similar, and 
quite similar to what was seen in the Sine-Gordon model; however when we looked at the 
details of the scattering we did spot some differences. The scattering in the first model 
was more elastic; in each case we related this to the fact that the solitons had different rest 
mass energies at the obstruction. The origin of such mass difference in the two models 
is due to the change of the soliton slope as the soliton climbs the obstruction. Due to 
this difference the scattering of the soliton on the obstruction was more elastic in model 
1 than in model 2 

When we looked at the scattring by the hole we found that, like in the Sine- Gordon 
case, the scattering was much more inelastic in both models with model 2 generating much 
more radiation. In fact, in model 2, for holes sufficiently deep the soliton was completely 
trapped in them; ie we could not find any velocity above which the soliton could escape 
from the hole. 

The question the arises whether all these results are very generic, ie hold in most 
(1+1) dimensional solitonic models or are just the property of the \ip 4 model. This 
problem is currently under investigation. Moreover, we need also to obtain some analytical 
understanding of the observed results. In the Sine-Gordon model this was achieved [2] by 
the introduction of effective variables describing the oscillation of the vacuum in the hole 
- very much in the way Goodman et al [5] explained the results on Fei et al [I]. We plan 
to look for a similar explanation in this case too. 

References 

[1] B.M.A.G. Piette, W.J. Zakrzewski and J. Brand, J. Phys. A 38, 10403-10412 (2005) 

[2] B.M.A.G. Piette and W.J. Zakrzewski, Scattering of Sine-Gordon kinks on potential 
wells, preprint to appear in J. Phys. A (2007) 

[3] Kurosh Javidan, J. Phys. A 39 10565-10574 (2006) 

[4] Z. Fei„ Yu.S. Kivshar, and L. Vazquez, Phys. Rev. A 45, 6019-6030 (1992) 
[5] R.H. Goodman and R. Haberman, Physica D 195, 303-323 (2004). 



11 



